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THE MILKY WAY HAS NO THICK DISK 

Jo Bow 1 ' 2,3 , Hans- Walter Rix 4 , & David W. Hogg 4 ' 5 

Astrophys. J., in press 

ABSTRACT 

Different stellar sub-populations of the Milky Way's stellar disk are known to have different vertical 
scale heights, their thickness increasing with age. Using SEGUE spectroscopic survey data, we have 
recently shown that mono-abundance sub-populations, defined in the [a/Fe]-[Fe/H] space, are well 
described by single exponential spatial-density profiles in both the radial and the vertical direction; 
therefore any star of a given abundance is clearly associated with a sub-population of scale height h z . 
Here, we work out how to determine the stellar surface-mass density contributions at the solar radius 
Ro of each such sub-population, accounting for the survey selection function, and for the fraction 
of the stellar population mass that is reflected in the spectroscopic target stars given populations of 
different abundances and their presumed age distributions. Taken together, this enables us to derive 
Sfl Q (/i z ), the surface-mass contributions of stellar populations with scale height h z . Surprisingly, we 
find no hint of a thin-thick disk bi-modality in this mass-weighted scale-height distribution, but a 
smoothly decreasing function, approximately £/? (/i 2 ) oc exp(— h z ), from h z w 200 pc to h z « 1 kpc. 
As h z is ultimately the structurally defining property of a thin or thick disk, this shows clearly that the 
Milky Way has a continuous and monotonic distribution of disk thicknesses: there is no 'thick disk' 
sensibly characterized as a distinct component. We discuss how our result is consistent with evidence 
for seeming bi-modality in purely geometric disk decompositions, or chemical abundances analyses. 
We constrain the total visible stellar surface-mass density at the Solar radius to be £ R = 30 ± 1 Mq 
pc . 

Subject headings: Galaxy: abundances — Galaxy: disk — Galaxy: evolution — Galaxy: formation 
— Galaxy: fundamental parameters — Galaxy: structure 



1. INTRODUCTION 

In both the Milky Way (Gilmore & Reid 1983) and 
external galaxies (Tsikoudi 1979; Burstein 1979), the 
stellar-density or luminosity profiles perpendicular to 
the disk plane reveal an excess over a simple expo- 
nential "thin disk" at distances \z\ > 1 kpc. This 
excess, confirmed by many studies in the last two 
decades (e.g., Reid & Majewski 1993; Juric et al. 2008), 
seems to be a generic feature of galactic disks (e.g., 
Yoachim & Dalcanton 2006) and has typically been de- 
scribed by, and often ascribed to, the presence of a 
"thick-disk" component, distinct from both the halo and 
thin disk components (e.g., Juric et al. 2008). However, 
whether or not the excess is sensibly attributed to a dis- 
tinct "thick-disk" component is, in our view, an open 
question (see § 1.1). 

Studies of the age, kinematics, and elemental- 
abundance ratios of probable members of the thicker 
disk component have revealed that it is older 
(e.g., Bensby et al. 2005), kinematically hotter (e.g., 
Chiba & Beers 2000), and metal-poor and enhanced in 
a elements (e.g., Fuhrmann 1998; Prochaska et al. 2000; 
Lee et al. 2011). Recent analyses (Navarro et al. 2011; 
Lee et al. 2011) also revealed a striking bi-modal dis- 
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tribution of disk stars in the elemental-abundance pa- 
rameter plane of [a/Fe] and [Fc/H]. However, these 
analyses did not account for proper volume sampling, 
and the [a/Fe] distribution will appear bi-modal even 
if the underlying (enrichment) age distribution is per- 
fectly smooth (Schonrich & Binney 2009a), simply be- 
cause [a/Fe] strongly changes as soon as enrichment by 
Type la supernovae (SNe la) becomes important. 

These conceptual issues can be illustrated by the re- 
cent finding that geometric decompositions of the Galac- 
tic disk yield strikingly different results for the struc- 
tural parameters (e.g., the radial scale lengths) for 
the thicker and thinner components, when they are 
based on an elemental-abundance selection of sample 
stars (Bensby et al. 2011; Bovy et al. 2011, Bll here- 
after), as opposed to luminosity- or volume-selection 
(e.g., Juric et al. 2008). This is because sample selection 
by elemental abundance, e.g., [a/Fe] and [Fe/H], can link 
sub-sample members independent of their position and 
velocity, while selections based only on luminosity, vol- 
ume or kinematics correlate the sample selection with the 
structure to be inferred. Therefore, discerning whether 
the thin- and thick-disk components should be thought of 
as distinct has been difficult on the basis of the existing 
evidence, given the modest number statistics and, most 
importantly, the conceptual or practical selection effects 
of many observational studies. Such issues have troubled 
the question of a thin-thick disk dichotomy for a long 
time (e.g., Nemec & Nemec 1991, 1993; Ryan & Norris 
1993; Norris 1999). 

The question of a thick-disk component plays an im- 
portant part in our understanding of how our Galaxy 
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formed and evolved. Various qualitatively dissimi- 
lar scenarios have been proposed to explain the pres- 
ence of a thick-disk component. Many of these in- 
voke mechanisms external to the already existing disk 
that may be expected in hierarchical structure forma- 
tion: direct accretion of stars from a disrupted satellite 
galaxy (Abadi et al. 2003), heating of a pre-existing thin 
disk through a minor merger (e.g., Quinn et al. 1993; 
Kazantzidis et al. 2008), or star formation induced by a 
gas-rich merger (e.g., Brook et al. 2004). However, quies- 
cent internal dynamical evolution can also reproduce the 
locally observed properties of the Milky Way's thick disk 
component (e.g., Schonrich & Binney 2009a,b). Recent 
attempts to differentiate these scenarios through volume- 
averaged stellar kinematics have not been conclusive 
(e.g., Sales et al. 2009; Dierickx et al. 2010; Wilson et al. 
2011). 

Discerning the mechanism that leads to the thick- 
ness distribution of galactic disks is important in con- 
straining the rate of minor mergers, the resilience of 
stellar disks to such mergers, and the importance of 
'internal' heating mechanisms such as radial migra- 
tion (Sellwood & Binney 2002) or turbulent, gravitation- 
ally unstable disk evolution (e.g., Bournaud et al. 2009; 
Forster Schreiber et al. 2009). While external heat- 
ing mechanisms or distinct disk-formation epochs may 
lead to sensibly distinct 'thin' and 'thick' disk com- 
ponents, internal mechanisms should not. Ultimately, 
these questions are also important for understanding 
how much of an archaeological record our Galactic disk 
holds and how much formation memory has been erased 
(Freeman & Bland-Hawthorne 2002). 

1.1. Approaches to dissecting the Galactic disk 

In our view there exists no prior analysis that can an- 
swer the question of whether it is sensible to view the 
thinner and thicker parts of the Milky Way's disk as dis- 
tinct components rather than a smooth distribution of 
stellar-disk scale heights. This is mainly for three rea- 
sons: 

• Either analyses have started out with insisting that 
two vertically-exponential components be fit to the ob- 
servations; this may be a sensible approach if the data 
are not good enough to warrant more complex mod- 
els, but it pre-supposes the answer. Bll showed that 
chemically defined sub-populations correspond to simple, 
single-exponential, sub-components of the Milky Way's 
disk. Bll also found a wide range of vertical scale height 
for these chemically-defined sub-populations, smoothly 
varying from thinner to thicker components with increas- 
ing (enrichment) age. 

• Decompositions have been based on geometric or 
kinematic sample definitions. Absent a very clear sep- 
aration of components in position or velocity space, this 
approach inevitably has an air of circular reasoning that 
precludes unique decompositions. Indeed, our recent 
analysis (Bll) has shown that state-of-the-art geomet- 
ric decompositions (e.g., Juric et al. 2008) cannot even 
reliably tell whether the thicker disk components are ra- 
dially more extended or more concentrated, as a compar- 
ison with a structure-independent abundance-based sub- 
sample selection has shown. Studies of the elemental- 
abundance trends of kincmatically-selected samples of 
stars (e.g., Prochaska et al. 2000; Feltzing et al. 2003) 



cannot assess the distinctness of the thick-disk compo- 
nent without correcting for the (strong) biases induced 
by the selection. The kinematics and spatial structure of 
a population of stars arc inextricably coupled through the 
dynamical properties of the population, and arc therefore 
not a priori independent. 

• Conceptually the right approach seems to be to define 
sub-samples by a property that may correlate with disk 
structure but is formally independent of it, such as stel- 
lar age or stellar abundances. Such 'tags' are properties 
that do not change along the orbit of the star, nor do they 
change if a star changes its orbit due to a minor merger or 
radial migration (Freeman & Bland-Hawthorne 2002); in 
this sense, tagging stellar populations by elemental abun- 
dances is even better than using the best structural or dy- 
namical tags available: integrals of motion, or orbital ac- 
tions. With such sub-samples, one can then ask whether 
their respective disk structures suggest distinct thin and 
thick disks. This approach has been re-advocated in re- 
cent studies (Navarro et al. 2011; Lee et al. 2011), draw- 
ing on the seemingly bi-modal distribution of stars in the 
[a/Fe] - [Fe/H] plane of elemental abundances; however, 
these studies did not account for selection effects and 
volume corrections. 

Believing that the last approach is the conceptually 
correct one, we derive in this paper for the first time a 
scale-height distribution of the stellar mass at the So- 
lar radius in the Galaxy to see whether this distribution 
shows any evidence of bi-modality. This requires two 
steps: first, one needs to associate a vertical scale-height 
with any given star. Our recent analysis (Bll) showed 
that sub-populations of stars of a given [a/Fe] and [Fe/H] 
are well described by a single exponential in their ver- 
tical profile, making ([a/Fe], [Fe/H]) a suitable tag for 
each star that uniquely determines the scale-height of its 
sub-population. Second, one needs to determine what 
the (surface) mass fraction of stars is in this [a/Fe]- 
[Fe/H] sub-population. As such samples can only be de- 
fined through spectroscopic surveys, this in turn requires 
proper accounting for the spectroscopic survey selection 
function (Bll), and an estimate of the stellar mass frac- 
tion that the spectroscopically targeted stars represent 
of the entire population at this elemental abundance. 

We work out this last step in this paper (§ 2), and 
on that basis derive the distribution of stellar surface- 
mass density- weighted distribution of stellar scale heights 
at i?o, ^R (h z ) (§ 2). Surprisingly, we find a smooth, 
exponentially-declining surface-mass-density spectrum 
as a function of scale height between « 200 and 1000 pc. 
This mass spectrum shows no gaps, excess, or hints of 
bi-modality at large scale heights, leading us to conclude 
that the Milky Way has no distinct thick disk. We dis- 
cuss implications, and explain how this relates to recent 
work that arrived at qualitatively different conclusions, 
in § 3. 

2. THE MASS- WEIGHTED SCALE HEIGHT DISTRIBUTION 
OF DISK STARS IN THE MILKY WAY 

Recently, we have performed (in Bll) number-density 
fits to sub-populations of stars defined as narrow boxes 
in the elemental-abundance plane spanned by metal- 
licity [Fc/H] and a-enhancement [a/Fe], based on G- 
dwarf spectra and star counts from SDSS / SEGUE 
(Abazajian et al. 2009; Yanny et al. 2009; SEGUE 
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Fig. 1. — Distribution of G-type dwarfs in the SDSS / SEGUE spectroscopic sample in the [Fe/H], [a/Fe] elemental-abundance plane. The 
left panel shows the raw number counts, while the right panel shows mass-weighted number counts that are corrected for spectroscopic 
selection effects and converted into total stellar surface-mass densities at the Solar radius using the stellar population modeling described 
in Appendix A. Pixels span 0.1 dcx in [Fe/H] and 0.05 dex in [a/Fe], which is larger than in the equivalent Figure in Lee et al. (2011), as 
the spatial number-density fitting used as part of the mass-weighting demands larger bins. The mass weighting shows that the prominent 
bi-modality seen in the raw number counts (left, cf. Lee et al. 2011) is mostly a consequence of the uneven spectroscopic sampling of the 
underlying stellar populations. The remaining hint of [a/Fe] bi-modality seen in the right panel is a natural consequence of the enrichment 
physics, i.e., SN la enrichment delays, even for a perfectly uniform (enrichment) age distribution (e.g., Schonrich & Binney 2009a). 



elemental-abundance uncertainties are 0.2 dex in [Fe/H] 
and 0.1 dex in [a/Fe]). G-type dwarfs are the most lumi- 
nous tracers whose main-sequence lifetime is larger than 
the expected disk age at basically all metallicities. We 
refer the reader to Bll for a discussion of the data set 
and for a detailed description and discussion of these 
fits, accounting for the volume correction of the spec- 
troscopic survey. Bll show that each mono-abundance 
sub-population has a simple density structure that can be 
described by an exponential profile both in Galactocen- 
tric radius R and vertical height z. The inferred scale 
heights for different mono-abundance sub-populations 
vary with [a/Fe] and [Fc/H], increasing smoothly from 
about 200 pc to about 1200 pc when going from popula- 
tions with near-solar abundances to populations that are 
more metal-poor and enhanced in a elements. 

These results were derived from the spectroscopic 
SEGUE G-dwarf sample, whose face- value distribution 
in the [Fe/H] -[a/Fe] plane is suggestive of two distinct 
populations definable in this elemental-abundance space, 
because there is a peak near solar abundances and one at 
metal-poor and a-enhanced abundances (see Lee et al. 
2011, Bll, and the left panel of FIG. 1). However, 
SEGUE targeted stars at high latitudes in an apparent- 
magnitude range 14.5 < r < 20.2 that finds G-type 
dwarfs at vertical heights > 500 pc, significantly above 
the bulk of the thin-disk population (which Bll inferred 
to have a scale height of 250 pc). Additionally, G-type 
dwarfs of different metallicities have different luminos- 
ity and hence survey volumes, and the spectroscopic 
targeting was weighted toward the fainter end of the 
apparent-magnitude range and toward higher latitudes 



(see Yanny et al. 2009 and the discussion in Appendix A 
of Bll). Thus, the distribution of the SEGUE sample 
in the [Fe/H]-[a/Fc] plane reflects the survey selection 
function more than the abundance distribution of the 
Milky Way stars within 1 to 2 kpc from the Sun. 

To properly assess the contribution of the various 
mono-abundance sub-populations to the mass or surface- 
mass budget in the Solar neighborhood, we must convert 
the number of spectroscopically-observed stars in each 
bin to a surface-mass density at the Solar neighborhood. 
On the one hand, this requires incorporation of a model 
for the SEGUE selection function (see Appendix A of 
Bll) and our exponential-disk fits (Bll). On the other 
hand, this requires the use of stellar-population models 
that can relate the observed number of G-type dwarfs 
to the total mass of the stellar population, given the 
metallicity of the sub-population and an assumed star 
formation history for it. This is described in detail in 
Appendix A. 

Briefly, for G-type dwarfs in a given ([Fe/H], [a/Fe]) 
bin we calculate their total number per square pc (in- 
tegrated over vertical height) by adjusting the normal- 
ization of the number-density profile such that after 
running it through our model for the SEGUE selec- 
tion function it predicts the observed number of stars 
in each bin. Then we relate the number density of 
G-type dwarfs to the total stellar surface-mass density 
by multiplying the number density by the average mass 
of a G-type dwarf — calculated using Padova isochroncs 
(Marigo et al. 2008) — and dividing it by the fraction of 
the mass in a stellar population in G-type dwarfs (calcu- 
lated using the same isochrones and assuming a lognor- 
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Fig. 2. — Distribution of stellar surface-mass density at the So- 
lar radius S^ (/i z ) as a function of vertical scale height h z . The 
thick black histogram shows the total stellar surface-mass density 
in bins in h z , calculated by summing the total stellar masses of sub- 
populations in bins in [a/Fe] and [Fe/H], The stellar surface- mass 
densities of the individual elemental-abundance bins in [Fe/H] and 
[a/Fe] are shown as dots, with values for Eij ([Fe/H], [a/Fe]) on 
the j/-axis. The points are color-coded by the value of [o/Fe] in 
each bin and the size of the points is proportional to the square 
root of the number of data points that the density fits are based 
on. Some of the errorbars are smaller than the points. Elemental 
abundance bins have a width of 0.1 in [Fe/H] and 0.05 in [ct/Fe]. 

mal Chabrier (2001) initial mass function). At a given 
abundance, this fraction of course depends on the age of 
the population, and here is calculated by marginalizing 
over a flat age distribution between 0.5 and 10 Gyr for 
each bin. However, averaging only over older ages for 
a-enhanced stars would be appropriate, as a-enhanced 
stars likely represent the oldest part of the disk. As we 
show in Appendix A, this gives similar results, with a 
slightly steeper decline in S(i?o) with h z . We calculate 
uncertainties on the surface-mass densities by varying the 
density-parameters according to the posterior probability 
distribution for the parameters in Bll. These uncertain- 
ties do not include systematic uncertainties due to the 
use of the stellar isochrones. This procedure results in 
an estimate of the stellar surface-mass density contribu- 
tion at the Solar radius for any abundance-selected sub- 
population, which in turn has a vertical scale height asso- 
ciated with it. The relative total stellar surface-mass den- 
sities of different mono-abundance bins are not affected 
by assuming a different initial mass function (IMF), al- 
though assuming a different IMF can systematically shift 
all surface-mass densities by a few percent (see below). 

The results from this mass estimation are shown in 
FIG.s 1 and 2. The left panel of FIG. 1 is simply a 
more coarsely binned version of the unweighted SEGUE 
G-dwarf sample abundance distribution, which shows 
two distinct maxima, one considerably more metal poor 
and a-enhanced than the other, seemingly reflecting a 
chemically-distinct thick-disk component. It is impor- 
tant to note that the marginalized [Fe/H] metallicity 
distribution (left panel, top) shows no hint of any bi- 
modality. There is distinct bi-modality in the marginal- 
ized [a/Fe] distribution, but Schonrich & Binney (2009a) 
already showed that even for a smooth age distribu- 



tion such bi-modality arises, simply separating stars 
that formed before and after enrichment by SN la be- 
came important. The right panel of FIG. 1 shows the 
stellar surface-mass density at the solar radius Rq in 
each elemental-abundance bin, corrected for selection ef- 
fects due to the spectroscopic SEGUE selection as de- 
scribed above. It represents the properly mass-weighted, 
underlying distribution of disk stars in the elemental- 
abundance space spanned by [Fe/H] and [a/Fe], and dra- 
matically differs from the raw sample distribution in the 
left panel: it does not have the strong bi-modality appar- 
ent in the raw SDSS / SEGUE number distribution. It 
also shows no hint of a bi- modal [Fc/H] metallicity dis- 
tribution and the remaining hint of [a/Fe] bi-modality is 
explained as in the left panel. 

The right panel of FIG. 1 now provides the rel- 
evant weights of each mono-abundance bin to a 
surface-mass-weighted distribution of disk scale heights. 
The colored symbols in FIG. 2 show exactly these 
surface-mass density contributions for all these mono- 
abundance bins in [Fe/H] and [a/Fe] (with widths 0.1 
in [Fe/H] and 0.05 in [a/Fe]) versus the scale height 
of those sub-population from Bll, color-coded by their 
[a/Fe] enhancement. We then sum the surface- mass 
contributions Efl (/i z |[Fc/H], [a/Fc]) of sub-populations 
([Fe/H], [a/Fe]) into bins in scale height. This results 
in the thick black histogram, which represents Y,R (h z ), 
or more simply p(h z ), the surface-mass weighted distri- 
bution of vertical scale heights in the Solar neighbor- 
hood (which we will refer to in the remainder simply as 
the 'scale-height distribution'). That is, for any random 
stellar-mass element this function gives the probability 
density for the scale height of the structural component 
to which it belongs. This is the function we set out to 
construct in order to examine whether it makes sense 
to think of distinct thin and thick disk components in 
the Milky Way. Remarkably, we find that the scale- 
height distribution simply decreases quite smoothly to- 
wards larger scale heights, with an approximately expo- 
nential relation between surface-mass density and scale 
height S(i?o) oc exp(— h z ). The scale height distribution 
does not show any gaps, excesses, or hints of bi-modality, 
beyond this simple relation. 

By combining all of the stellar surface-mass density 
estimates we can precisely measure the total visible stel- 
lar surface-mass density at the Solar radius. We find 
S flo = 30 ± 1 Mq pc -2 . This is similar to the estimate 
of Flynn et al. (2006), who report S^ o = 29 M pc" 2 . 
This estimate depends slightly on the assummcd IMF. 
Using the exponential IMF (IMF3) of Chabrier (2001) 
gives T,* Ro = 29.5 M Q pc -2 ; the IMF from Chabrier 
(2003) gives Y,* Ro = 29 M s pc" 2 ; and a Kroupa (2003) 
IMF gives Z Ro = 32 M pc" 2 . 



3. DISCUSSION 

FIG. 1 shows that properly correcting for the spectro- 
scopic sampling of the underlying stellar sub-populations 
is crucial in assessing the elemental-abundance distribu- 
tion at i?o of spectroscopically selected samples of stars. 
The abundance distribution without this correction is 
heavily influenced by the survey-specific spatial and mass 
sampling of the underlying stellar population — both of 
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which act to make the metal-poor and a-enhanced sub- 
populations more prominent in the high-latitude and 
color-selected SEGUE sample — which leads to a spu- 
riously enhanced bi-modality in elemental-abundance 
space. The mass-weighted metallicity distribution in 
FIG. 1 has no bi-modality. The mass- weighted [a/Fe]- 
distribution in the same Figure has only a hint of a bi- 
modality, as is expected in standard smooth star for- 
mation and enrichment scenarios (Schonrich & Binney 
2009a), and reflects enrichment physics and not galaxy 
evolution. 

Improper weighting of the distribution of structural pa- 
rameters can again lead to spurious bi-modal signatures. 
FIG. 5 in Bll shows the location of chemically-defined 
sub-populations in the space of the structural param- 
eters (radial scale length, vertical scale height). This 
scatterplot, based on equal-area bins in ([Fc/H],[a/Fc]), 
gives undue prominence to the low-[Fc/H], high-[a/Fc] 
(with respect to solar abundances) bins, as many of them 
only contribute a negligible amount to the total stellar 
mass. Here, FIG. 2 shows that the proper mass- weighting 
of the ([Fe/H],[o!/Fe]) sub-populations gives a vertical- 
scale-height distribution that is smooth and monoton- 
ically declining between thinner disk component scale 
heights of 200 pc and thicker components' scale heights of 
w 1200 pc, with no gaps or excesses beyond the smooth, 
approximately exponential distribution. Bll found that 
each elemental-abundance bin was preferentially fit by a 
single exponential rather than by two disk components, 
such that the smooth scale-height distribution in FIG. 2 
is not merely the result of the smoothing out of an in- 
trinsically bi-modal distribution by elemental-abundance 
errors (which are small for the SEGUE sample, see § 2) 
or overlapping abundance distributions of distinct "thin" 
and "thick" disks; if either of these were the case Bll 
should have detected two components in the abundance 
bins with single-exponential scale heights in the range of 
approximately 400 to 600 pc. The large uncertainties on 
the radial scale lengths of the chemically-defined mono- 
abundance sub-populations in Bll complicate a similar 
assessment of the radial structure of the disk, but this 
has no bearing on the analysis of the vertical structure. 
Upcoming surveys such as APOGEE (Eisenstein et al. 
2012) or HERMES/GALAH (Freeman et al. 2010) that 
will sample stellar populations in the plane of the Milky 
Way will be able to study the radial structure in more 
detail. 

Thus, stars in the Solar neighborhood have a smoothly 
decreasing probability of belonging to structural compo- 
nents with increasing scale heights. This implies that the 
thicker disk component in the Milky Way is simply the 
tail of a continuous and monotonic scale-height distribu- 
tion. This has been suggested before (e.g., Norris 1987; 
Schonrich & Binney 2009b) , but never directly measured 
as we do here. As such, there is no distinct thick-disk 
component in our Galaxy. 

Together with the findings in Bll that the thicker and 
older components of the Galactic disk have a shorter ra- 
dial scale lengths than the thinner and younger compo- 
nents this qualitatively points toward a continuous inter- 
nal mechanism such as radial migration or turbulent disk 
evolution being predominantly responsible for the thick- 



ening of the disk, rather than an external merger or heat- 
ing event. However, a rigorous comparison with models, 
where thick stellar-disk components arise from one or a 
few distinct events triggered, e.g., by satellite infall, is 
needed to see whether the present data are indeed incon- 
sistent with the results presented here. Formation age 
of stars may serve as a sensible marker in simulations to 
associate individual stars with a 'parent sub-population' 
whose scale height can be determined. In this context, 
it is worth noting that combined data from APOGEE 
and Kepler will soon provide ages for a large number 
of stars through asteroseismology (e.g., Gilliland et al. 
2010), which could help in mapping mono-age popula- 
tions into the mono-abundance populations studied here. 

The proper selection function analysis and subsequent 
conversion to stellar mass densities can of course be used 
in the future to look with proper mass weighting at 
other disk diagnostics beyond the spatial distribution, 
such as the orbital eccentricity or vertical motions, which 
has not been done correctly (e.g., Dierickx et al. 2010; 
Wilson et al. 2011). 

Our results that show that the Milky Way has no dis- 
tinct thick-disk component might appear to be at odds 
with observations of external edge-on galaxies, where 
"thick disk" components are found to be universal (e.g., 
Yoachim & Dalcanton 2006). However, those decompo- 
sitions are restricted to luminosity-weighted geometric 
decompositions and suffer from the uncertain influence of 
dust in the mid-plane of the galaxies; they cannot isolate 
components based on elemental abundances; and they 
only perform discrete two-component fits, which natu- 
rally prefers two components over a single exponential 
component, even though the underlying scale-height dis- 
tribution might be more complicated, as shown in this 
paper. The findings from external galaxies might be 
more correctly interpreted as proving the universality of 
thicker disk components rather than that of "thick disks" . 
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APPENDIX 

THE SURFACE-MASS DENSITY IN A MONO- ABUNDANCE SUB-POPULATION 

In practice, we calculate the surface-mass density contribution at the Solar radius from each mono-abundance sub- 
component, defined as a box in the [Pe/H]— [a/Fe] plane, by first converting the observed number counts into a 'number 
column density' (N(Rq) in the notation of Bll; units stars pc -2 ) and then converting this number density into a 
total stellar surface-mass density using stellar-population models. This involves estimating the mean (individual) 
stellar masses of the stars we have in the sample, and estimating which fraction of their total stellar population at 
that abundance that mass range constitutes. 

When performing the stellar number-density fits in Bll we did not fit for the normalization of the density. However, 
we can calculate this normalization by adjusting it such that it predicts the observed number of stars when we run 
the density model through the observational selection function. Thus, we calculate the predicted number of stars for a 
density normalized to have unit surface-mass density (in stars pc -2 ); the correct normalization constant is then the 
total number of observed stars divided by the predicted number of stars for unit surface-mass density. We calculate the 
uncertainty in this number by performing this procedure for each of the samples from the posterior distribution for the 
number-density parameters, obtained in Bll to estimate the uncertainties in the best-fit number-density parameters. 
This gives for each abundance bin a number N(R ; [Fe/H], [a/Fe]) ± cat ; [f c /h],[q/Fo] that represents the total number 
of stars per square pc in this [Fe/H] -[a/Fe] bin at the Solar radius. 

In order to turn this number density N(Rq; [Fc/H], [a/Fe]) into a stellar mass we use stellar isochrones in the SDSS 
photometric system (Girardi et al. 2004; Marigo et al. 2008; Girardi et al. 2010) 6 to model the mass M{g — r, r, [Fc/H]) 
as a function of color g — r, age r, and metallicity [Fc/H]. The isochrones all assume [a/Fe] = 0. A lognormal Chabrier 
(2001) mass function 4>(M) is used to weight the contribution of various masses to the stellar population's mass (i.e., 
the mass function provides the measure in mass space). Assuming different forms for the IMF systematically changes 
the stellar surface-mass densities by a few percent independent of metallicity. 

The total stellar mass is related to the number of G-type dwarfs approximately as 

stellar mass in a [Fc/H] - [a/Fc] bin = N(R ; [Fe/H], [a/Fe]) (M G )([Fc/H]) w _1 ([Fc/H]) , (Al) 

where (Afc}([Fe/H]) is the average mass of a G-type dwarf in the 0.45 < g — r < 0.58 color range 7 and u;([Fe/H]) 
is the ratio of a stellar population's mass in G-type dwarfs to the total mass in the population. We calculate both 
(Mg}([Fc/H]) and w([Fe/H]) using the stellar isochrones and assuming a lognormal Chabrier (2001) initial mass 
function. 

We calculate the mass in G-type dwarfs of a stellar population as 

M G ,totai([Fc/H]) = f d(g - r) [ "** dr f dM cf>(M) M(g - r, r, [Fc/H]) , (A2) 

JO. 45 Jr mm J 

where we thus use a flat prior in age and we only use the dwarf part of the isochronc. The total stellar mass of a stellar 
population is given by a similar expression, but without the color restriction. In what follows we use (r m i n ,r max ) = 
(0.5,10) Gyr for all abundance bins. As [a/Fe] is a relative age indicator (see discussion in Bll) with populations with 
[a/Fe] > 0.25 probably > 7 Gyr old an [a/Fe]-dcpcndent age distribution would be more appropriate, although the 
exact form it should take is hard to establish. We also used (r min , r max ) = (1,8) Gyr for [a/Fe] < 0.25 and (r min , r max ) 
= (7,10) Gyr for [a/Fe] > 025 and found no significant difference in the inferred stellar masses for all of the abundance 
bins. w([Fe/H]) is then the ratio of the stellar mass in G-type dwarfs to the total stellar mass. The ratio w([Fe/H]) is 
approximately given by 

w([Fe/H]) « 0.0425 + 0.0198 [Fe/H] + 0.0057 [Fe/H] 2 , (A3) 

between -1.5 < [Fc/H] < 0.3. 

We calculate the average mass of a G-typc dwarf in our color range using a similar procedure. We calculate the 
average mass by marginalizing over the same age distribution, but without weighting by the mass function as the 
color dependence of the mass function is very limited over this narrow color range and the distribution of the SEGUE 
G-dwarf sample is uniform in color g — r. The average mass is thus calculated as 

<A/ G )([Fe/H]) = f 58 d(g - r) f dr M(g - r, r, [Fc/H]) ; (A4) 

JO. 45 Jr min 

this is approximately given by 

(M G /M )([Fe/H]) w 0.956 + 0.205 [Fe/H] + 0.051 [Fc/H] 2 , (A5) 



again between -1.5 < [Fe/H] < 0.3. 

6 Retrieved using the Web interface provided by 
Leo Girardi at the Astronomical Observatory of Padua 
http : //stev . oapd. inaf . it/cgi-bin/cmd_2 . 3 

7 Here and in equations (A2) and (A4) we use a slightly wider 



color range than that of the SDSS / SEGUE G-dwarf sample, as 
the Padova isochrones do not have very high resolution in g — r 
color. 
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The average mass (Mc)([Fe/H]) of a G-type dwarf as a function of [Fc/H] only changes by about 20 percent when 
going from metal-poor to metal-rich stars, with more metal-poor stars having smaller masses. The fractional contri- 
bution to the mass budget of G-type dwarfs w([Fc/H]) decreases from 5 percent for metal-rich stars to approximately 
3 percent for metal-poor stars. Thus, these mass-correction factors are not the main drivers for the structure in FIG. 2. 
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